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ABSTRACT 

Context. Solar-like oscillations have now been observed in several stars, thanks to ground-based spectroscopic observations and 
space-borne photometry. CoRoT, which has been in orbit since December 2006, has observed the star HD49933 twice. The oscillation 
spectrum of this star has proven difficult to interpret. 

Aims. Thanks to a new timeseries provided by CoRoT, we aim to provide a robust description of the oscillations in HD49933, i.e., to 
identify the degrees of the observed modes, and to measure mode frequencies, widths, amplitudes and the average rotational splitting. 
Methods. Several methods were used to model the Fourier spectrum: Maximum Likelihood Estimators and Bayesian analysis using 
Markov Chain Monte-Carlo techniques. 

Results. The different methods yield consistent result, and allow us to make a robust identification of the modes and to extract precise 
mode parameters. Only the rotational splitting remains difficult to estimate precisely, but is clearly relatively large (several fiHz in 
size). 
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1. Introduction 

Stars are now the objects of seismic studies after decades of 
similar studies for the Sun, thanks to the advent of space-borne 
photometric observations (e.g., MOST, CoRoT and Kepler) and 
extremely precise ground-based spectroscopic ob servations (for 
a complete review, see for e.g. lAerts et al.ll2008l) . This applies 
in particular to stars presenting solar-like p modes (acoustic 
oscillations stochastically excited by convection), with CoRoT 
observations of such st ars showing clearly individual peaks in 
the Fourier spectra (e.g. Mic hel et al.ll2008l) . Among these stars, 
HD49933 has already been the target of asteroseismic cam- 
paigns. It was first observed spec troscopically from the ground 
for 10 nights dMosser et al.l|2005l) . In 2007, a first photometric 
time series of 60 days was collected by CoRoT, followed by a 
new long run of 137 days in 2008. HD49933 is a F5 main se- 
quence star with an apparent visual magnitude my = 5.77. It 



* The CoRoT space mission, launched on 2006 December 27, was 
developed and is operated by the CNES, with participation of the 
Science Programs of ESA, ESA's RSSD, Austria, Belgium, Brazil, 
Germany and Spain. 



is hotter than the Sun (r eff =6780K or r eff =6 500 K, iBruntt et all 
120081: lBmndl2009r. iRvabchikova et al Il2009h with an estimated 
mass of ~ 1.2M (Mosser et al. 2005) and an estimated radius 
of 1.34 + 0.06 fl Q dThevenin et alJ 120061) . The su rface rotation 
(vsinQ was determined to be around lOkms 1 dMosser et al.l 
l2005tlSolano" et al.l [2005h . The surface rotation period has also 
been measured at ~ 3.4 days, using the 60-day Co RoT timeseries 
dAppourchaux et aL1l2008t iDeheuvels e t al. 2008), from the sig- 
natures of photospheric transiting active regions (e.g., spots) 
which give rise to a clear peak in the very low-frequency part 
of the Fourier spectrum. 

The seis mic interpreta t ion of HD49933 has proven to be 
very difficult. Mos ser et alJ ([2005) could not isolate individual p 
modes in the Fourier spectrum of observed line-of-sight veloci- 
ties but were able to find a regular pattern in the spectrum, which 
is the signature of the large frequency separation between modes 
of same degree / (but increasin g radial order n). The first CoRoT 
time-series was analysed by Appourchau x et alJ d2008l) . and 
these data clearly show individual p-mode peaks in the Fourier 
spectrum. However, the peaks show large widths, making the 
interpretation less than straightforward: a given peak could be 
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interpreted as being a closely spaced pair of /=0 and 1—2 modes, 
or a single (but rotationally split) /=1 mode. Based on the mod- 
eling of the sp ectrum using a Max i mum L ikelihood Estimator 
fitting method, Appourchau x et a Q d2008l) chose one of these 
two possible interpretations (hereafter called model A) based 
on the highest likelihood of eac h model. This first time-s eries 
was the object of other studi es. | Appourchaux et al.| C 2009) put 
into perspective the results of lAppourchaux et alj d2008l) . show- 
ing that the likelihood ratio test does not give the probability of 
the hypothesis given the data, but only the sig nificance of the 
data given the hypothesis. iBenomar et alj (|2009), who applied a 
Bayesian analysis to the same time series, could not definitely 
favour one interpretation (model A) over the alternate (model 
B), based on the whole probability distribution of each model. 
Gruberbauer etalJ (2009), using a Ba yesian approach too , also 
consider the identification ambiguous. lGaulme et al.1 (|2009) used 
a simpler Bayesian approach (Maximum A Posteriori, or MAP, 
approach). The most probab le model they found corre sponds 
to the s ame identification as lAppourchaux et alj d2008l) . More 
recently, Mosser & Appourchaux ( 2009J) proposed an empirical 
method to determine the identification of the modes. Its appli- 
cation considers the model B as the more likely when using the 
two datasets used here. 

It should be noted that the case of HD49933 is quite different 
from the solar case: solar modes are very narrow in comparison. 
Their widths (1 //Hz for the modes with the highest amplitudes) 
are much smaller than the small-frequency separations between 
1=0 modes of order n and the neighbouring 1=2 modes of order 
n-\ (being typically around 10//Hz for the Sun). Moreover, the 
star inclination angle, if small, tends to attenuate the visibility of 
mode components with azimutal order m + 0, making the mode 
identification even more difficult. 

Here, we use the new long-run CoRoT observations of 
HD49933, together with the original shorter run timeseries 
(Fig-EK in order to properly describe the acoustic oscillations 
of the star clearly visible in the Fourier spectrum (Fig.Q] Fig. [5] 
and Fig.©. 

2. Methodology 

T he time series used here w ere extracted in the same way as 
in lAppourchaux et al. (2008). The gaps in the series represent 
slightly less than 10% and were filled by linear interpolation. In 
order to provide a robust seismic interpretation of the CoRoT 
timeseries for HD49933, several analyses were performed, us- 
ing different methods (Maximum Likelihood Estimators, or 
MLE; and various Bayesian analyses); or the same method ap- 
plied in an independent manner. These methods were already 
used in the same way to analyze the initial run of 60 days 
(lAppourchaux et al. 2008t IBenomar et alj|2009l) . The "fitters" of 
Appourchaux et al. (2008) and some additional fitters sought to 
find a best-fitting model spectrum. The model has several con- 
tributions. There is a contribution from the background, which 
includes signatures of convection and possibly phenomena with 
longer time scales (e.g. those related to the stellar activity) and 
contributions from the individual p modes. Each mode is de- 
scribed by a set of parameters: a central frequency, a width and a 
height. A single height and width was fitted to each I = 0/2 pair 
and the closest I = 1 mode (in frequency). The relative heights 
of the 1=0, 1 and 2 modes took fixed values, which were as- 
sumed to be independent of frequency. A single rotational fre- 
quency splitting parameter was fitted to all non-radial modes. 
The stellar inclination, which governs the relative heights of the 
different m components for a given (n,l) mode, was also fitted 



as a single, global parameter. The observed spectrum used by 
the different fitters was an averaged spectrum (frequency reso- 
lution 0.1 9 //Hz) made from three timeseries of the same dura- 
tion, which came from the two CoRoT runs: the first 60-day run, 
IRaOl, and two 60-day long timeseries from the longer second 
run, LRaOl (data available at idoc-CoRoT.ias.u-psud.fr). 
The different fitters analysed this spectrum independently, and 
we then compared the results. 

3. Results 

The availability of the new longer CoRoT timeseries makes the 
mode (degree identification significantly less ambiguous than 
it was before. The p-mode peaks are still observed to be very 
wide (several //Hz), and the typical signal-to-noise ratio (SNR) 
(defined as the ratio of the height of a mode to the level of the 
background around the mode) is similar to that for the first, 
shorter CoRoT run. However, the additional information pro- 
vided by the longer second run is sufficient to allow the modes 
to be tagged with far greater confidence than was hitherto pos- 
sible. There was a very good agreement between the results of 
the fitters, with model B (which has an /=1 mode at 1755 //Hz; 
see Table[T]for details) favoured strongly over Model A. As part 
of the analysis for this paper, w e performed a Bayesian analysis 
(using the method described by IBenom ar 2008; Benoma r et alj 
2009) to compare model A with model B, and this now favours 
model B at a confidence level higher than 99 %. It is important to 
add that this confidence level does depend to some extent upon 
the hypothesis used for the modeling (e.g. on the a priori con- 
straints applied to the model parameters). Here, the 99 % level 
was estimated for a model with fixed mode height ratios (see 
Section|2] above). When the fixed height ratio constraint was re- 
laxed so that individual heights were fitted to modes of differ- 
ent I, we found that the confidence level for model B dropped 
no lower than 95 %. (It is worth adding that the average fitted 
height ratio agrees, to within errors, with the solar-like value we 
adopted as the fixed height ratio.) 

We also established, using the Bayesian approach, that mod- 
els that include I = 2 modes are strongly favoured over those 
which do not (at a confidence level over 99.9 %). To check for 
evidence in the Fourier spectrum of I = 3 modes, we instead used 
a collapsed spectrum. We could find no evidence for a significant 
excess of power in the wings of the / = 1 modes. Evidently, the 
SNR in the / = 3 modes is too low for them to be observed. 

We make one final remark r egarding the identification prob- 
lem. lAppourchaux et ail (120081) . who favoured model A, made 
their choice based solely on a comparison of the maximum like- 
lihoods given by a classical MLE analysis. When the fitting 
problem is non-trivial, this type of analysis can converge on a 
subsidiary maximum (not the true, global maximum), biasing 
any statistical comparison of two possible models. In addition 
to MLE, we also applied for the work in this paper the Markov 
Chain Monte-Carlo (MCMC) analysis, which circumvents the 
above problem by giving a complete sampling of the parame- 
ter space. Both MCMC and MLE now clearly favour Model B. 
The present results are due to the conjunction of two facts. First, 
the extra information provided by the new observations (adding 
137 days to the first 60 days of observation). Second, the con- 
firmation given by the use of Bayesian analysis about a robust 
comparison of the associated probabilities to each model. All 
this ensures the identification without ambiguity, whatever the 
method used. 

Most of the frequencies for model B returned by the dif- 
ferent fitters lie within a l-cr error interval (which corresponds 
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Fig. 1. Mean power spectrum using 3 time series of 60 days (solid grey line), and the fitted model (dashed line). 
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Fig. 2. Large frequency separation computed from 1=0 (crosses 
linked by a dashed line) and 1=1 (diamonds linked by a solid 
line) modes. 



to a level of confidence of 68%), and all do so within a 1.5- 
cr interval. These error intervals are particularly small in the 
frequency range v € [1300, 2000] /iHz where the modes have 
relatively large heights. Agreement is particularly good at 1=1 
(cr ~ 0.6/iHz; see Fig. 9). This is because these modes are not af- 
fected by prominent nearby modes, as is the case for the closely 
spaced / = and 2 modes (which show significant overlap in 
frequency). The I = 2 modes have lower amplitudes than their 
I = neighbours and consequently have the largest errors of any 
of the observed / (<x ~ 2yuHz, see Fig. 9). 

Table Q] lists Model B frequencies from one of the analy- 
ses used here (a Ba yesian analysis cou pled to MCMC sampling; 
see for example iBenomar et al.l l2009). The echelle diagram of 
these frequencies is shown in Fig. [3] From the most reliable 1=0 
and 1=1 modes, it is possible to estimate the large frequency 
separations Ao and Ai and their variation with frequency (see 
Fig. [2). The uncovered frequency variation may be regarded as 
being significant, given the good precision on the estimated fre- 
quencies. Any estimate of the frequency difference 602 between 



neighbouring 1=0 and 1=2 modes is much less reliable because 
of the difficulty of fitting these modes (see above). Thus, the only 
result that can be provided is the average: (£02) = 4.7/iHz. 

The widths of the modes, which are related to the damping, 
are also listed in Table [T] and shown in Fig. [4] This figure illus- 
trates the main difficulty of interpreting the spectrum: the rela- 
tively high values of the widths. As on the Sun, an increase in 
width with increasing frequency is visible, but the level of pre- 
cision prevents a more detailed analysis. The amplitudes of the 
modes, which depend on the balance between damping and ex- 
citation, are listed in Table Q] The / = amplitudes are plotted 
in Fig. [6] The maximum is about ~3.7 ppm, which, while higher 
than the maximum amplitude of low-/ solar oscillations, is still 
lower than predicted by a scalin g of amplitude on {L/M) a (with 
a =s 0.7, see Samadi et al l2007h . Extraction of the rotational fre- 
quency splitting, v s , of the modes remains very difficult: the anal- 
yses of the different fitters failed to converge on a unique solu- 
tion. Estimated values were in the range 3.5^/Hz < v s < 6.0^Hz, 
i.e., high compared to the solar value, but as expected given the 
surface rotation period estimated from the low-frequency part of 
the Fourier spectrum (P mt ^3.4 days). Finally, we were able to 
extract robust values for the angle of inclination of the star. All 
analyses converged on an angle of of 17° +7. This is in agree- 
ment with an independent determination m ade using measure- 
ments of the stellar v sin 2, radius and period (ISolano et alj 2005; 
Mos ser et al.l20 09). However, as mentioned earlier, such a small 
angle greatly favours the visibility of m = mode components, 
rendering the splitting measurement very difficult. 

4. Conclusion 

The two C0R0T observation runs on the star HD49933 - the first 
60-day run, and the more recent longer run - have now provided 
enough data to resolve the identification of modes in the oscilla- 
tion spectrum at a very high confidence level. This identification 
relied on several independent analyses. The new data have also 
allowed us to improve the precision in the mode parameters, with 
fractional improvements being in the range from 40 % to 70 % 
depending on the parameter. It is now possible to determine pre- 
cise mode frequencies for /=0 and 1 modes. The 1=2 mode pa- 
rameters are more difficult to estimate because of the overlap 
with the stronger, neighbouring 1=0 modes. The widths and am- 
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Fig. 3. Echelle diagram built from the fitted frequencies for a 
large separation of 85 //Hz. The error boxes indicate the lcr (68% 
confidence) interval. 
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Fig. 4. Widths of the fitted modes, with lcr error bars 



plitudes of the modes are well determined, as is the inclination 
angle of the star. The rotational frequency splitting remains the 
only mode parameter that is poorly constrained. However, it is 
clearly much higher than the solar value, and similar or larger 
in size to the inverse of the surface rotation period. In summary, 
this new information on the acoustic oscillations of HD49933 
now opens the possibility for detailed seismic modeling of the 
star. 

Acknowledgements. W.J.C. and Y.E. wish to thank the UK Science and 
Technology Facilities Council (STFC) for support under grant ST/F00204/1. 
I.W.R. and G.A.V. also thank STFC, for support under grant PP/E001793/1. J.B. 
acknowledges support through the ANR project Siroco. We thank J. Leibacher 
for very useful comments. 



References 

Aerts, C, Christensen-Dalsgaard, J., Cunha, M., & Kurtz, D. W. 2008, 
Sol. Phys.,251, 3 

Appourchaux, T, Michel, E., Auvergne, M, et al. 2008, A&A, 488, 705 
Appourchaux, T., Samadi, R., & Dupret, M. A. 2009, A&A, in press 
Benomar, O. 2008, Communications in Asteroseismology, 157, 98 
Benomar, O., Appourchaux, T., & Baudin, F. 2009, A&A, in press 
Bruntt, H. 2009, A&A, accepted 

Bruntt, H., De Cat, P., & Aerts, C. 2008, A&A, 478, 487 



Degree 


Frequency 


Itr 


O 


Height/Noise 


/ 


(pHz) 


interval 


interval 


ratio 





1206.25 


+2.23/-5.57 


1 no / "7 1~ 

+3.92/-7 .26 


0.6 





1 o o o rv7 

1288.97 


.1 1 C / f\ f\f\ 

+ 1.25/-0.90 


. on/ 1 ryi 

+3.20/-1.97 


0.7 





1 O Tl 1 A 

1373.34 


+0.93/-1.16 


1 ha / cn 

+ 1.70/-2.50 


1.3 





1460.14 


+0.65/-0.80 


1 OO /OIC 

+ 1.22/-2.15 


1.4 





1544.69 


i n o c / n nc 

+0.85/-0.95 


+ 1. /6/-Z.0J 


1.2 





1631.10 


+0.59/-0.69 


+ 1.17/- 1.59 


2.4 





1714.49 


+ 1.06/-1.17 


1 nn / 1 
+2.00/-2.21 


1.9 





1799.75 


i t\ O A 1 1 AT 

+0.84/- 1.07 


1 /"O / O /"O 

+ 1.62/-2.62 


2.4 





1884.82 


+0.88/-1.26 


+ 1.63/-2.65 


2.4 





1972.73 


i C\ /ATI 

+O.0//-O./I 


i 1 on/ 1 c n 

+ 1.30/-1.59 


2.3 





2057.82 


. O OO / 1 O /I 

+3.22/-1.34 


A ^70 / O O A 

+4.72/-2.24 


1.9 


u 


Zl^f / . 1U 


_i_ n 7 A/ H 
+U. /0/-U.OJ 


-1-1 A^/ 1 QQ 
+ 1 .OJ/-1 .yy 


1 A 


U 


2ZJo.4o 


+Z. 1 1/-Z.83 


1 O O C / C yl A 

+ 3.oj/-j.44 




u.y 


U 


Ti oo i n 
z3zZ.ll) 


i o ti / o m 
+2.22/-3.03 


1 A no / C CA 

+4.08/- J .04 


U. / 


U 


o<inc ca 
2408. JO 


+ 1.9//-2.33 


1 1 H C / /f OI 

+3. 0/-4.81 


n c 
U.5 


U 


O/IOC 7A 
Z4V J . / O 


+3.34/-3.21 


i 7 nn/ a 00. 
+ / .00/-0. 23 


U.3 


U 


lean oc 


i /i in/ i in 
+4. /0/-3.2V 


i t £i / 1 n 7 
+ /.0I/-IO. / 


U.J 


1 


1 161.54 


. n o c / a m 

+0.85/-0.91 


. 1 ^ "7 / O O C" 

+ 1.67/-2.25 


0.6 


1 


1244.63 


i 1 no / 1 n 
+ 1. 02/- 1.1 / 


1 on/ ^ c 


0.7 




1 TOO 1/1 

1328.34 


+0.70/-0.65 


1 1 /in/ 1 £L 
+ 1.40/-1.26 


1.2 


1 


1414.93 


. n c"7 / n co 

+0.57/-0.58 


1 00 / 1 00 
+ 1.22/- 1.28 


1.4 


1 


1500.54 


. n HA / A TO 

+0.70/-0.78 


+ 1.36/-1.62 


1.2 


1 


1586.62 


+0.48/-0.49 


1 n no / 1 no 

+0.98/- 1.02 


2.3 


1 


1670.48 


. n c"7 / n c o 

+0.57/-0.58 


+ 1.1 6/- 1.23 


1.9 


1 


1755.30 


. n c i / n co 

+0.51/-0.53 


1 t\/~ 1 1 f\ c 

+ 1.06/- 1.05 


2.3 


1 


1840.68 


i n a n / n cn 

+0.49/-0.50 


1 n no / 1 no 

+0.98/- 1.02 


2.3 


1 


1 no o 1 1 

1928.13 


i n cn / n c 1 

+0.50/-0.51 


1 1 nn / 1 n c 

+ 1.00/- 1.05 


2.3 


J 


2014.38 


i n c /i / n c /i 

+0.54/-0.54 


1 1 m / 1 nn 

+ 1.07/- 1.09 


1.8 




o i ni Co 
Z1U1.JO 


, n £7/ n 70 
+U.0 //-U. /Z 


,1 /io/ 1 CA 
+ 1 .4Z/-1 . JO 


1 A 


1 


1 1 nn o 1 
21V0.81 


1 n nn/ n nn 

+o.yo/-o.90 


1 1 00/1 no 
+ 1.88/-1.92 


0.8 


1 


11T7 on 

22/ /. 89 


+ 1.1 6/- 1.14 


1 O O O / O T7 

+2.33/-2.37 


0.7 


1 


Tito 

2362. /o 


i 1 A 1 / 1 A 1 
+ I.OI/-I.OI 


+J.39/-3.48 


U.5 




o /i en q c 
24j0.3j 


+2.20/-2. / 1 


+4.D8/-/.01 


U.J 


1 


iron 

2j39.4V 


1 1 cn/ A C 
+ 1.j0/-4.3j 


+ 3.24/-1 1.2 


U.3 


2 


1 199.91 


1 1 no / c n /i 
+4.08/- 5 .04 


1 / 1 1 1 

+8. / //-1 1.3 


0.6 


2 


1 O "7 A 

1287.24 


.O CI/ 1 T7 

+3.51/-3.77 


. "1 CO/ O /" c 

+7.53/-8.65 


0.7 


2 


1369.60 


1 O AO 1 1 11 

+2.48/-3.1 1 


1 C O C / H 7H 

+5.25/- 7 .79 


1.3 


2 


1455.42 


1 1 11 / 1 C\A 

+2.33/-1.94 


1 c /i n / a r\H 

+5.49/-4.07 


1.4 


2 


1541.54 


+3.0//-4.4/ 


1 c n a 1 n /i 

+5.94/-9.48 


1.2 


2 


1626.30 


. in / i aa 

+2.79/-3.00 


+5.16/-6.17 


2.4 


2 


1712.67 


1 O 1 / O A A 


1 C O O / /I H C 

+5.23/-4.75 


1.9 


2 


1794.39 


1 71 / 1 n 
+2. /3/-2. 10 


1 /n / n 1 

+6.47/-3.91 


2.4 


2 


1 O O 1 o o 

1881.83 


. ha / no 

+2.70/-2.02 


. c ot / -1 no 

+5.27/-4.08 


2.4 


2 


1965.19 


1 n 1 / 1 h a 
+2.01/-1. /4 


1 a no 1 c c 

+4.78/-3.55 


2.3 


2 


2060.22 


. O O O / C 1 ~7 

+2.88/-5.17 


A CI 1 1 An 

+4.57/-7.40 


1.9 


7 




_i_r! 7c/ 9 f.7 


_i_7 A?/ ^ AS 

/ XlLj- J .uo 


1 A 


2 


2230.68 


+5.10/-2.91 


+9.99/-5.89 


0.9 


2 


2316.86 


+4.35/-3.35 


+9.06/-8.24 


0.7 


2 


2403.52 


+3.97/-4.56 


+7.77/- 11.1 


0.5 


2 


2491.50 


+4.96/-5.16 


+10.3/- 11.2 


0.3 


2 


2576.51 


+5.03/-7.94 


+9.92/- 18.6 


0.3 



Table 1. Frequencies of the fitted modes. The 1 and 2-cr intervals 
correspond to confidence levels of 68% and 95%, respectively. 
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Fig. 5. Echelle diagram computed for a large separation of 85 /iHz for the mean power spectrum. The diagram is smoother to 0.8 
//Hz. The identification of the ridges is far from evident. The most probable model correspond to I — 1 for the left ridge and I = 
for the right ridge (same convention as in Fig. [3]) 
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Table 2. Amplitudes and widths of the fitted modes. The 1-cr intervals correspond to confidence levels of 68%. 
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Fig. 7. Lightcurves of HD49933 for the 60-day long initial run (left) and the 137-day long run (right) 




1620 1640 1660 

Frequency (/iHz) 



1 680 



Fig. 8. Detail of the spectrum (solid grey line), and of the fitted model (dashed line) with a pair I = 0/2 on the left and an / = 1 on 
ther right. 
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Fig. 9. Mean (upper and lower) lcr error bars for the frequencies of 1=0 (triangles), 1=1 (diamonds) and 1=2 (crosses) modes, 
showing the decreasing precision on the frequencies from 1=1 to 1=2 modes. 



